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We present a simultaneous measurement of the forward-backward asymmetry and the top-quark 
polarization in tt production in dilepton final states using 9.7 fb _1 of proton-antiproton collisions 
at y/s = 1.96 TeV with the DO detector. To reconstruct the distributions of kinematic observables 
we employ a matrix element technique that calculates the likelihood of the possible tt kinematic 
configurations. After accounting for the presence of background events and for calibration effects, 
we obtain a forward-backward asymmetry of A u = (15.0±6.4 (stat) ±4.9 (syst))% and a top-quark 
polarization times spin analyzing power in the beam basis of kP = (7.2 ± 10.5 (stat) ±4.2 (syst))%, 
with a correlation of —56% between the measurements. If we constrain the forward-backward 
asymmetry to its expected standard model value, we obtain a measurement of the top polarization 
of reP = (11.3±9.1 (stat)±1.9 (syst))%. If we constrain the top polarization to its expected standard 
model value, we measure a forward-backward asymmetry of A u = (17.5 ± 5.6 (stat) ± 3.1 (syst))%. 
A combination with the DO A u measurement in the lepton±jets final state yields an asymmetry of 
A u = (11.8 ± 2.5 (stat) ± 1.3 (syst))%. Within their respective uncertainties, all these results are 
consistent with the standard model expectations. 

PACS numbers: 14.65.Ha, 12.38.Qk, 13.85.Qk, 11.30.Er 


I. INTRODUCTION 

In proton-antiproton collisions at i/s=1.96 TeV, top 
quark pairs are predominantly produced in valence 
quark-antiquark annihilations. The standard model 
(SM) predicts this process to be slightly forward- 
backward asymmetric: the top quark (antiquark) tends 
to be emitted in the same direction as the incoming quark 
(antiquark), and thus, in the same direction as the incom¬ 
ing proton (antiproton). The forward-backward asym¬ 
metry in the production is mainly due to positive con¬ 
tributions from the interference between tree-level and 
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next-to-leading-order (NLO) box diagrams. It receives 
smaller negative contributions from the interference be¬ 
tween initial and final state radiation. The interferences 
with electroweak processes increase the asymmetry. In 
the SM, the asymmetry is predicted to be ~ 10% (U- 
H|. Within the SM, the longitudinal polarizations of the 
top quark and antiquark are due to parity violating elec¬ 
troweak contributions to the production process. The 
polarization is expected to be < 0.5% for all choices of 
the spin quantization axis am. 

Physics beyond the SM could affect the tt production 
mechanism and thus both the forward-backward asym¬ 
metry and the top quark and antiquark polarizations. In 
particular, models with a new parity violating interaction 
such as models with axigluons [a-Q) can induce a large 
positive or negative asymmetry together with a sizable 
polarization. 

The tt production asymmetry, A u , is defined in terms 
of the difference between the rapidities of the top and 
antitop quarks, A y ti = y t - yy. 

A ti = N(Ay tl > 0) - N(Ay ti < 0) 

N(Ay t t > 0) ± N(Ay t t < 0) ’ 

where N(X) is the number of events in configuration X. 
By definition, A u is independent of effects from the top 
quark decay such as top quark polarization. However, it 
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requires the reconstruction of the tt initial state from the 
decay products, which is challenging especially in dilep¬ 
ton channels. 

Measurements of A u have been performed in the lep- 
ton+jets channels by the CDF H and DO [III Collabo¬ 
rations. Other asymmetry measurements have been per¬ 
formed using observables based on the pseudo-rapidity 
of the leptons from t —> Wb —> £ub decays fl2l4l5l ]. All 
these measurements agree with the SM predictions. A 
comprehensive review of asymmetry measurements per¬ 
formed at the Tevatron can be found in Ref. H- 

As top quarks decay before they hadronize, their spin 
properties are transferred to the decay products. The 
top (antitop) polarization Pf (Pr) along a given quan¬ 
tization axis n impacts the angular distribution of the 
positively (negatively) charged lepton [1] 

A^ = I(l + * ± P,W), (2) 

where 9 + (9~) is the angle between the positively (neg¬ 
atively) charged lepton in the top (antitop) rest frame 
and the quantization axis n, and k + (k~) is the spin 
analyzing power of the positively (negatively) charged 
lepton, which is close to 1 (—1) at the 0.1% level within 
the SM |||. The polarization terms n + P/[ ( k ~P//) can 
be obtained as two times the asymmetry of the cos 9 + 
(cos# - ) distribution 

N(cos9 ± > 0) — N(cos9 ± < 0) 

" A"(cos 9 ± > 0) + N (cos 9 ± < 0) 

In the following we use the beam basis, where n is the 
direction of the proton beam in the tt zero momentum 
frame. Since we only use the beam basis, we omit the 
subscript n in the following and define the polarization 
observable as: 

kP = i (k+P + - k~P~) = A e+ - A r . (4) 

Polarization effects have been studied at the Tevatron 
in the context of the measurements of the leptonic asym¬ 
metries in Ref. H , but no actual measurement of the 
polarization has been performed. Measurements of the 
polarization have been conducted for top pair produc¬ 
tion in pp collisions at the Large Hadron Collider at 
y/s = 7 TeV. These measurements, performed in differ¬ 
ent basis choices, are all consistent with the SM expec¬ 
tations EH. 

This article presents a simultaneous measurement of 
A tt and kP with the DO detector in the dilepton de¬ 
cay channel. It is based on the full Tevatron integrated 
luminosity of 9.7 fb -1 using tt final states with two lep¬ 
tons, ee, e/z, or pp. We first reconstruct the A y t i and 
cos distributions employing a matrix element integra¬ 
tion technique similar to that used for the top-quark 
mass measurement in the dilepton channel [20]. These 
distributions are used to extract raw measurements of 


asymmetry and polarization, A[ aw and «:P raw , in data. 
The experimental observables A[ aw and reP raw are corre¬ 
lated because of acceptance and resolution effects. Us¬ 
ing a MC@NLO Ej 0] simulation, we compute the re¬ 
lation between the raw measurements AH,„ and 
and the true parton-level asymmetry and polarization to 
determine calibration corrections. We then extract the 
final measured values of A u and kP. This is the first 
measurement of the tt forward-backward asymmetry ob¬ 
tained from the reconstructed Adistribution in the 
dilepton channel and the first measurement of the top 
quark polarization at the Fermilab Tevatron collider. 


II. DETECTOR AND OBJECT 
RECONSTRUCTION 

The DO detector used for the Run II of the Fermi¬ 
lab Tevatron collider is described in detail in Refs. [23l — 
(H- The innermost part of the detector is composed 
of a central tracking system with a silicon microstrip 
tracker (SMT) and a central fiber tracker embedded 
within a 2 T solenoidal magnet. The tracking system 
is surrounded by a central preshower detector and a 
liquid-argon/uranium calorimeter with electromagnetic, 
fine hadronic, and coarse hadronic sections. The central 
calorimeter (CC) covers pseudorapidities [27] of \rj\ < 
1.1. Two end calorimeters (EC) extend the coverage to 
|? 7 | < 4.2, while the coverage of the pseudorapidity re¬ 
gion 1.1 < \rj\ < 1.5, where the EC and CC overlap, is 
augmented with scintillating tiles. A muon spectrome¬ 
ter, with pseudorapidity coverage of |? 7 | < 2, is located 
outside the calorimetry and comprises drift tubes and 
scintillation counters, before and after iron toroidal mag¬ 
nets. Trigger decisions are based on information from the 
tracking detectors, calorimeters, and muon spectrometer. 

Electrons are reconstructed as isolated clusters in the 
electromagnetic calorimeter and required to spatially 
match a track in the central trac king system. They have 
to pass a boosted decision tree El criterion based on 
calorimeter shower shape observables, calorimeter iso¬ 
lation, a spatial track match probability estimate, and 
the ratio of the electron cluster energy to track momen¬ 
tum (E/p). Electrons are required to be in the accep¬ 
tance of the electromagnetic calorimeter (|t?| < 1.1 or 
1.5 < \r]\ < 2.5). 

Muons are identified by the presence of at least one 
track segment reconstructed in the acceptance (|7y| < 2.0) 
of the muon spectrometer that is spatially consistent with 
a track in the central tracking detector Hi- The trans¬ 
verse momentum and charge are measured by the curva¬ 
ture in the central tracking system. The angular distance 
to the nearest jet, the momenta of charged particles in a 
cone around the muon track, and the energy deposited 
around the muon trajectory in the calorimeter, are used 
to select isolated muons. 

Jets are reconstructed from energy deposits in the 
calorimeter using an iterative midpoint cone algo- 
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rithm (30] with a cone radius of 72. = 0.5 [31,]. The jet 
energies are calibrated using transverse momentum bal¬ 
ance in 7 +jet events [32|. 


III. DATASET AND EVENT SELECTION 

The signature of tt production in dilepton final states 
consists of two high-p^ leptons (electrons or muons), two 
high-pr jets arising from the showering of two b quarks, 
and missing transverse energy (Ifir) due to the undetected 
neutrinos. The main backgrounds in this final state arise 
from Z —> ££, with £ = e, n, or r, and diboson pro¬ 
duction (WW, WZ , ZZ ). These backgrounds are eval¬ 
uated from Monte Carlo (MC) simulated samples as de¬ 
scribed in section II V Cl Another source of background 
comes from fV+jets and multijet events, if one or two 
jets are misreconstructed as electrons or if a muon from 
a jet passes the isolation criteria. The contribution from 
these backgrounds, denoted as “instrumental background 
events”, are estimated directly from data as described in 
section HV El Each of the dilepton channels is subject to a 
different mixture and level of background contamination, 
in particular for the background arising from the Z -A ll 
process. We therefore apply slightly different selection 
requirements. The main selection criteria to obtain the 
final samples of tt candidate events are: 

1. We select two high px (pt >15 GeV) isolated lep¬ 
tons of opposite charge. 

2. We require that at least one electron passes a sin¬ 
gle electron trigger condition in the ee channel 
(r; 100 % efficient), and that at least one muon 
passes a single muon trigger condition in the pp 
channel (« 85% efficient). In the ep channel, we 
do not require any specific trigger condition, i.e., 
we use all DO trigger terms (« 100% efficient). 

3. We require two or more jets of px > 20 GeV and 
M < 2.5. 

4. We further improve the purity of the selection by 
exploiting the significant imbalance of transverse 
energy due to undetected neutrinos and by exploit¬ 
ing several topological variables: 

(i) The missing transverse energy ]fix is the mag¬ 
nitude of the missing transverse momentum, 
obtained from the vector sum of the trans¬ 
verse components of energy deposits in the 
calorimeter, corrected for the differences in de¬ 
tector response of the reconstructed muons, 
electrons, and jets. 

(ii) The missing transverse energy significance, 

is the logarithm of the probability to 
measure Ifix under the hypothesis that the 
true missing transverse momentum is zero, 


accounting for the energy resolution of indi¬ 
vidual reconstructed objects and underlying 
event 133; . 

(iii) Ht is the scalar sum of transverse momenta 
of the leading lepton and the two leading jets. 

In the ee channel we require > 5, in the ep 
channel Ht > 110 GeV, and in the pp channel 
> 5 and $t > 40 GeV. 

5. We require that at least one of the two leading jets 
be 5-tagged, using a cut on the multivariate dis¬ 
criminant described in Ref. [34[ . The requirement is 
optimized separately for each channel. The tt selec¬ 
tion efficiencies for these requirements are ~ 82%, 
« 83%, and ss 75% for the ee, ep,, and pp channels, 
respectively. 

6. The integration of the matrix elements by VEGAS, 
described in section IV Al may return a tiny proba¬ 
bility if the event is not consistent with the tt event 
hypothesis due to numerical instabilities in the in¬ 
tegration process. After removing low probability 
events, we retain signal events in the MC simula¬ 
tion with an efficiency of 99.97%. For background 
MC, the efficiency is > 99.3%. We remove no data 
events with this requirement. 


IV. SIGNAL AND BACKGROUND SAMPLES 
A. Signal 

To simulate the tt signal, we employ MC events gen¬ 
erated with the CTEQ6M1 parton distribution func¬ 
tions (PDFs) [35]] and MC@NLO 3.4 [2l], [2i| interfaced 
to herwig 6.510 [36] for showering and hadronization. 
Alternate signal MC samples are generated to study sys¬ 
tematic uncertainties and the shape of the A y t j distribu¬ 
tion. We use a sample generated with ALPGEN (37] inter¬ 
faced to PYTHIA 6.4 [231 for showering and hadroniza¬ 
tion and a sample generated with ALPGEN interfaced 
to HERWIG 6.510. For both samples we use CTEQ6l1 
PDFs [H. 

The MC@NLO generator is used for the nominal signal 
sample as it simulates NLO effects yielding non-zero A u . 
The value of A u at parton level without applying any se¬ 
lection requirement is A tt = (5.23 ±0.07 (stat))%, which 
is smaller than a SM prediction @] that includes higher 
order effects. 

The MC events are generated with a top-quark mass of 
mt = 172.5 GeV. They are normalized to a tt production 
cross section of 7.45 pb, which corresponds to the calcu¬ 
lation of Ref. |39| for m t = 172.5 GeV. The generated 
top mass of 172.5 GeV differs from the Tevatron average 
mass ofl73.18±0.94 GeV [4(|. We correct for this small 
difference in section fVI Bl 
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B. Beyond standard model benchmarks 

We also study the five benchmark axigluons models 
proposed in Ref. 41] that modify tt production. For 
each of the proposed beyond standard model (BSM) 
benchmarks, we produce a tt MC sample using the 
MADGRAPH [4i| generator interfaced to PYTHIA 6.4 for 
showering and hadronization, and the CTEQ6l1 PDFs. 
The Z' boson model proposed in Ref. [JO] is not consid¬ 
ered here since it is excluded by our tt differential cross- 
section measurement [43 ]. 


In the e/z channel, the template for instrumental back¬ 
ground events is obtained with the same selection criteria 
as for the samples of tt candidate events, but without ap¬ 
plying the complete set of electron selection criteria. For 
the /z/z and ee channels, the contributions from instru¬ 
mental background events is negligible and the result is 
not sensitive to the choice of template. For simplicity, 
we re-employ the e/z template for both the /z/z and ee 
channels. 


F. Comparison of MC simulation to selected data 


C. Background estimated with simulated events 

The background samples are generated using the 
CTEQ6L1 PDFs. The Z —> i£ events are generated 
using ALPGEN interfaced to PYTHIA 6.4. We normalize 
the Z —> ££ sample to the NNLO cross section [44|- The 
Pt distribution of Z bosons is weighted to match the dis¬ 
tribution observed in data [4f|, taking into account its 
dependence on the number of reconstructed jets. The 
diboson backgrounds are simulated using PYTHIA and 
are normalized to the NLO cross section calculation per¬ 
formed by MCFM [ihl. l47j. 


D. DO simulation 

The signal and background processes except in¬ 
strumental background are simulated with a detailed 
GEANT3-based j48| MC simulation of the DO detector. 
They are processed with the same reconstruction soft¬ 
ware as used for data. In order to model the effects of 
multiple pp interactions, the MC events are overlaid with 
events from random pp collisions with the same luminos¬ 
ity distribution as data. The jet energy calibration is 
adjusted in simulated events to match the one measured 
in data. Corrections for residual differences between data 
and simulation are applied to electrons, muons, and jets 
for both identification efficiencies and energy resolutions. 


E. Instrumental background estimated with data 

The normalization of events with jets misidentified as 
electrons is estimated using the “matrix method” [49j 
separately for the ee and e/z channels. The contribution 
from jets producing identified muons in the /z/z channel 
is obtained using the same selection criteria as for the 
sample of tt candidate events, but demanding that the 
leptons have the same charge. In the e/z channel, it is 
obtained in the same way but after subtracting the con¬ 
tribution from events with jets misidentified as electrons. 

Once the absolute contribution of instrumental back¬ 
ground events has been determined, we also need “tem¬ 
plate samples” that model their kinematic properties. 


A comparison between the expected and observed 
numbers of events at the final selection levels is reported 
in Table ffl The selected sample is relatively pure with 
a background fraction varying between 10% and 16% 
depending on the channel. A comparison of kinematic 
distributions between data and expectations at the final 
selection level is shown in Fig. [Q 

V. MATRIX ELEMENT METHOD 

To reconstruct distributions of kinematic observables 
describing the tt events, we use a novel modification of 
the matrix element (ME) integration developed for the 
rrit measurements |20j, |50| by the DO Collaboration. In 
particular, this method is employed to reconstruct the 
A y t i, cos(0 + ), and cos(0 - ) distributions, from which an 
estimate of the forward-backward asymmetry and top po¬ 
larization are extracted. 


A. Matrix element integration 


The ME integration used in Refs. (uni consists in 
computing the likelihood L z to observe a given event with 
the vector of measured quantities z, 


L z 


= - X! / /pdf(9i)/pdf(<?2)IT {pt) 

A ■ o'tot „ . J 

flavors x ,q-± , Q2 i 

x W{x,z) —= h== d<& 6 dp t Td4> t *dq 1 dq2. (5) 
4^(91 • ®) 2 


In this expression, x is a vector describing the kinematic 
quantities of the six particles of the pp —I tt —I l + vb£~vb 
final state, ^ is the matrix element describing the dy¬ 
namics of the process, dd> 6 is the 6-body phase space 
term, the functions /pdf are the PDFs of the incoming 
partons of momenta q\ and </2 and of different possible 
flavors, W(x,z), referred to as the transfer function, de¬ 
scribes the probability density of a parton state x to be 
reconstructed as z, W(p t f) is a function describing the 
distribution of the tt system tranverse momentum, 
while the azimuthal angle of this system, i j> tt , is assumed 
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TABLE I: Comparison between expected and observed numbers of events at the final selection level for the different channels. 
The values are reported with their statistical uncertainties. 


Channel Z il Dibosons Instrumental tt lljj Total expected Data 


yy 

10 . 65+8 

i 7+u.i 

- 0.1 

O-O-oio 

7Q Q+U-b 

' y ‘°- 0.6 

917 +uy 

i _ 0 7 

92 

1.00 

± 

0.10 

e/i 

13.03/gJ 

q 7 + 0.2 
' - 0.2 

i6.4±g : ; 

283.1+g 

316 . 21 + 

346 

1.09 

± 

0.05 

ee 

12.92±°-| 

i Q+ 0 - 1 
± -^-o.i 

1 q + 0.08 
0.08 

95 . 5+8 

112 l +0 - 8 

0 8 

104 

0.92 

± 

0.10 


to have a uniform distribution over [0, 2n\, and A ■ <Tt 0 t is 
the product of the experimental acceptance and the pro¬ 
duction cross-section. The matrix element, . is com¬ 
puted at leading order (LO) for qq annihilation only, as 
it represents the main subprocess (« 85%) of the to¬ 
tal tt production. The functions /pdf are given by the 
CTEQ6L1 leading order PDF set. The function W(p + 
is derived from parton-level simulated events generated 
with ALPGEN interfaced to pythia. More details on this 
function can be found in Ref. [5lj| . Ambiguities between 
partons and reconstructed particle assignments are prop¬ 
erly treated by defining an effective transfer function that 
sums over all the different assignments As we consider 
only the two leading jets in the integration process, there 
are only two possibilities to assign a given jet to either 
the b or 6 partons. 

The number of variables to integrate is given by the six 
three-vectors of final state partons (of known mass), the 
tt transverse momentum and transverse direction, and 
the longitudinal momenta of the two incoming partons. 
These 22 integration variables are reduced by the fol¬ 
lowing constraints: the lepton and 6-quark directions are 
assumed to be perfectly measured (8 constraints), the 
energy-momentum between the initial state and the fi¬ 
nal state is conserved (4 constraints), the £ + v and 
system have a mass of M\y = 80.4 GeV j52j] (2 con¬ 
straints), and the i + vb and £~Db system have a mass of 
mt = 172.5 GeV (2 constraints). Transfer functions ac¬ 
count for muon and jet energies. The transfer functions 
are the same as used in Ref. [50]. The electron momen¬ 
tum measurement has a precision of ss 3%, which is much 
better than the muon momentum resolution of typically 
10% and the jet momentum resolution of typically 20%. 
We thus consider that the electron momenta are perfectly 
measured. This gives one additional constraint in the ey 
channel and two additional constraints in the ee channel. 
Thus, we integrate over 4, 5, and 6 variables in the ee, ey, 
and yy channels, respectively. The integration variables 
are <j) tt , energy of leading jet, energy of sub-leading 
jet, and energy of the muon(s) (if applicable). 

The integration is performed usin g th e MC-based nu¬ 
merical integration program VEGAS [5j, |54[ . The inter¬ 
face to the VEGAS integration algorithm is provided by 
the GNU Scientific Library (GSL) [55|. The MC inte¬ 
gration consists of randomly sampling the space of in¬ 
tegration variables, computing a weight for each of the 
random points that accounts for both the integrand and 
the elementary volume of the sampling space, and finally 
summing all of the weights. The random sampling is 


based on a grid in the space of integration that is itera¬ 
tively optimized to ensure fine sampling in regions with 
large variations of the integrand. For each of the random 
points, equations are solved to transform these integra¬ 
tion variables into the parton-level variables of Eq. ©, 
accounting for the measured quantities z. The Jacobian 
of the transformation is also computed to ensure proper 
weighting of the sampling space elementary volume. 

B. Likelihood of a parton-level observable 

For any kinematic quantity K reconstructed from the 
parton momenta x, for example K(x) ~ yt~ Vt, we can 
build a probability density L Z (K) that measures the like¬ 
lihood of K{x) at the partonic level to give the recon¬ 
structed value K. This likelihood is obtained by insert¬ 
ing a term 5(K(x) — K) in the integrand of Eq. ©, and 
normalizing the function so that f L z (K)dK = 1. The 
probability density is obtained by modifying the VEGAS 
integration algorithm. For each reconstructed tt event 
and each point in the integration space tested by VEGAS, 
the integrand of Eq. © and the quantity K are com¬ 
puted. After the full space of integration has been sam¬ 
pled, we obtain a weighted distribution of the variable 
K that represents the function L Z (K) up to an overall 
normalization factor. 

For each reconstructed event with observed kinematics 
Zi, where i is an event index, we obtain a likelihood func¬ 
tion L Zi {K). By accumulating these likelihood functions 
over the sample of events, we obtain a distribution that 
estimates the true distribution of the variable K. The 
performance of this method of reconstruction for parton- 
level distributions is estimated by comparing the accu¬ 
mulation of likelihood functions to the true parton-level 
quantities for MC events, as shown in Fig. [2] 

C. Raw estimate of A tf 

We could choose to use the maximum of the likelihood 
function L z (Ay t t) to estimate the true value of A y t t on 
an event-by-event basis. However, to maximize the use 
of available information, we keep the full shape of the 
L z functions and accumulate these functions over the 
sample of tt events to obtain an estimate of the par- 
ton level distributions, which is then used to determine 
A u . This method has been verified to perform better 
than the maximum likelihood method. The distribution 
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FIG. 1: [color online] Comparison of distributions between data and MC simulations at the final selection for (a) the transverse 
momentum of the leading lepton, (b) the transverse momentum of the secondary lepton, (c) the pseudorapidity of the leading 
lepton, (d) the pseudorapidity of the secondary lepton, (e) the transverse momentum of the leading jet, (f) the transverse 
momentum of the secondary jet, (g) the Ht, and (h) the difference between the two lepton pseudorapidities. The overflow bin 
content has been added to the last bin. 







































9 



FIG. 2: [color online] Accumulation of likelihood functions (^2 eventa L Zi (K), with K along the vertical axis) versus the corre¬ 
sponding true parton level quantity (ATtrue along the horizontal axis) in tt MC events after applying the selection criteria for (a) 
K = (b) K = cos 8 + , and (c) K = cos 8 . Each single MC event i contributes in these plots with a complete distribution, 

L Zi (X), along the vertical axis for a given value on the horizontal axis, K trae . The shades of color indicate the bin contents in 
arbitrary units. 


Eevents L *i ( ^Utt) is shown in Fig. (Ha), after subtract¬ 
ing the background contributions from the data. The 
raw asymmetry A‘ aw , extracted from this distribution, is 
reported in Table El Since this A y t t distribution is an 
approximate estimation of the true distribution of A y t t, 
the raw asymmetry A*. aw is an approximation of the true 
A tt . The measurement therefore needs to be calibrated. 
The calibration is discussed below. 


TABLE II: Raw forward-backward asymmetry in data be¬ 
fore background subtraction, A^ a , asymmetry of the back¬ 
ground, , and measurement once the background contri¬ 
bution has been subtracted, A*aw- Asymmetries are reported 
in percent, together with their statistical uncertainties. 


Channel 

A data 
^raw 

A bkg 

raw 

Att _ A data—bkg 

^raw ^raw 

e/i 

9.2± 3.8 

0.3± 1.9 

10.1 ±4.2 

ee 

15.8 ±6.4 

0.1 ± 2.0 

18.8 ±7.6 

pp 

6.7 ± 7.9 

—0.3± 3.3 

7.8± 9.1 

Dilepton 

10.1 ± 3.0 

0.1 ± 1.1 

11.3 ± 3.4 


The use of an event-by-event likelihood function allows 
us to define an asymmetry observable for each event 

A= L z (Ay t i)dAy t f — / L z (Ay ti )dAy th (6) 
J 0 J — oo 

where the observable A averaged over the sample of tt 
candidate events is equal to the raw asymmetry A‘ aw . 
By construction, A lies in the interval [—1,+1]. For a 
perfectly reconstructed event without resolution effects, 
A would be either equal to —1 for A y t j < 0 or to +1 
for A y t i > 0. The use of A allows us to determine the 
statistical uncertainty on A* aw as the uncertainty on the 
average of a distribution. 


D. Raw estimate of kP 

In the same way as in the previous section, we 
use the accumulation of the likelihoods L z (cos9 + ) 
and L z (cos9~) to estimate the distributions of cos0 + 
and cosThe distributions ^events L Zi (cos9 + ) and 
^events L Zi (cos9~) are shown in Figs.^b) andEJc), af¬ 
ter subtracting the background contributions from the 
data. The raw asymmetries, A[: aw and A* aw , and the 
raw polarization nP rav , = Af aw — A^ aw extracted from 
the data are reported in Table Hill As for Aj aw , the mea¬ 
surement of KP raw needs to be calibrated to retrieve the 
parton-level values of the polarization. 


E. Statistical correlation between Alt w and fcP raw 

We measure the statistical correlation between A* aw 
and K-P raw in the data, which is needed to determine the 
statistical correlation between the measurements of A u 
and kP. In the same way as A* aw is the average of an 
event-by-event asymmetry A, the raw asymmetries A^ aw 
and A^ aw are the averages of event-by-event asymmetries 
denoted by A cos g+ and A cos g-. The correlation between 
At* and is identical to the correlation between the 

observables A and (A cos g+ — A cos g~). This correlation 
is determined from the background subtracted data by 
computing the RMS and mean values of the distributions 
of A, (A cos g+ A cosg -) : and A • (A cos $+ ^cos#~)- 

cor (A‘* aw ,«;P raw ) = 

< A ■ (A cos g+ — A COS 0- ) > — A raw ■ AtP raw 

RMS(A) • RMS(A cos6) + — A cose -) ' [) 

We report the values measured in data in Table hyi 
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TABLE III: Asymmetry estimates for the cos 0 ± distributions. The raw asymmetry measurement in the data before background 
subtraction, Af av ; data , the asymmetry of the background, Af a ^ bkg , and the measurement once the background contribution has 
been subtracted, Af av ; data_bkg , are reported. The polarization estimates defined as kP^ = Af^, xx — Af av ; xx are also given. 
All values are reported in percent, together with their statistical uncertainties. 


Channel 

a ^, data 

■^raw 

Al + , bkg 
^raw 

, data—bkg 
■^raw 

a£ , data 
•^raw 

A e ,bkg 

7-1 raw 

, data—bkg 
^raw 

pdata 

raw 

pbkg 
rx,± raw 

K p _ K pd a t a —bkg 

r ' J± raw rv - t raw 

e/i 

5.7 ± 4.1 

0.6 ± 2.1 

6.2 ±4.6 

—3.3±4.1 

2.6 ±2.1 

-4.0 ±4.6 

9.0± 5.8 

-2.0 ±2.4 

10.2 ±6.4 

ee 

13.4 ±7.2 

-3.2 ±2.0 

16.5 ±8.6 

-0.8 ±7.2 

—0.5 ± 2.1 

-0.9 ±8.6 

14.2 ±10.1 

—2.7 ± 2.3 

17.4 ±12.0 

gg 

—9.4± 8.1 

3.9 ±3.6 

— 11.5 ± 9.4 

—3.7 ± 8.1 

2.3 ±3.5 

—4.7 ±9.3 

—5.7 ± 11.8 

1.5 ± 3.7 

-6.9 ±13.7 

Dilepton 

4.6 ±3.3 

0.2 ± 1.3 

5.2 ±3.7 

—2.9±3.3 

1.7 ± 1.2 

—3.5 ±3.7 

7.5 ±4.7 

— 1.5 ± 1.4 

8.7 ± 5.3 
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FIG. 3: [color online] Estimated distribution of the (a) At/ t f, (b) cos 9 + , and (c) cos# - observables in dilepton events 
after subtracting the expected background contribution. Deviations beweeen the background-subtracted data and MC can 
be attributed to statistical fluctuations. The background-subtracted data asymmetries and the MC asymmetries extracted 
from these distributions are also reported. These raw asymmetries need to be corrected for calibration effects to retrieve the 
parton-level asymmetries. 


TABLE IV: Measurement of the statistical correlation be¬ 
tween the asymmetry A( aw and the polarization /tP raw for the 
data, background, and background subtracted data. Values 
are reported in percent, together with their statistical uncer¬ 
tainties. 


Channel 

Data Background Data—Background 

e/x 

27 ± 6 

9± 3 

28 ± 6 

ee 

10± 12 

9± 3 

9 ±14 

gg 

36 ±10 

6± 5 

39 ±12 

Dilepton 

26 ± 5 

9± 2 

28 ± 5 


VI. RESULTS CORRECTED FOR 
CALIBRATION 

The calibration procedure finds a relation between the 
raw asymmetry and polarization, (A( aw , KP raw ), obtained 
after subtracting the background contributions, and the 
true asymmetry and polarization ( A tt 1 nP ) of tt events. 
The calibration procedure corrects for dilution effects 
that arise from the limited acceptance for tt events, the 
finite resolution of the kinematic reconstruction, and the 
simplified assumptions used in the matrix element inte¬ 
gration ( e.g ., leading order ME, no gg tt ME, only 
two jets considered). The relation is inverted to extract 


a measurement of A tl and nP from the values of A* aw 
and ftP raw observed in data. 

The nominal calibration is determined using a sample 
of simulated tt MC@NLO dilepton events. The procedure 
is repeated with the samples from the other generators 
(see section IIV Al and IIVBI) to determine different sys¬ 
tematic uncertainties. We normalize the individual ee, 
eg, and gg contributions to have the same proportions 
as observed in the data samples after subtracting the ex¬ 
pected backgrounds. 


A. Samples for calibration 


We produce test samples from a nominal MC sample 
by reweighting the events according to the true value of 
the parton-level A y t i, cos 6 + , and cosf? - . The reweight¬ 
ing factors are computed as follows. 
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1. Reweighting of lepton angular distributions 


The general expression for the double differential lep¬ 
ton angle distribution is [EJ 


d 2 a 

dcos6 + dcos6~ 


— (l + k + P + cos 9 + + k P cos 9 
CcosO + cos9~) , (8) 


where C is the spin correlation coefficient, which is 
~ 90% in the SM. In the beam basis one has nP ~ 
K+P+ PS —tv P . We use this relation to reweight a 
given MC sample to simulate a target polarization of 
K-Ptest = 3 ( K+P+ -K-P-). 


2. Reweighting of Ay t j distribution 

To determine a method of reweighting the A y t i distri¬ 
bution, denoted T>{Ay t i), we study its shape using the 
different ti MC samples of section IIVI at the generated 
level, i.e., before event selection and reconstruction. In¬ 
spired by the studies performed for the distribution of 
rapidity of the charged leptons in Ref. [56} . we rewrite 
2 ?(A y t i) as 


£*(A dti) = \ (D(Ay t i) + V(-Ay t t)) ■ (1 + A{Ay ti )), 

( 9 ) 

where A{Ay t i) = p(Ay“-)+-p(-Ay“-) is the ratio between 
the odd and even part of the Ay t j distribution, also called 
differential asymmetry as a function of Ay t j; we then fit 
A (Ay t t) with an empirical odd function 

/(Ayti) = P x tanh + P x , (10) 


where a and 7 are shape parameters, while /3 is a magni¬ 
tude parameter. The term /3 x was not needed 

in the study of Ref. [56| , but improves the modeling sig¬ 
nificantly for the case of A y t j. The results of the fit 
for different ti MC samples are shown in Fig. [I] If we 
reweight a MC sample so that the even part of the Ay t i 
distribution, the term a , and the term 7 are preserved, 
then the forward-backward asymmetry is proportional to 

P- 

These considerations yield the following procedure to 
produce a sample of test asymmetry starting from a 
MC sample of generated asymmetry %“mpie- We first fit 
the differential asymmetry at the generated level A(Ay t t) 
with the function / of Eq. m and determine the param¬ 
eters a, /?, and 7 . Then we apply weights to the events 
processed through the DO simulation 


w(A y tl ) 


1 + A t test /(Ay tt -) 

''sample 

1 + /(A y t i) 


( 11 ) 



FIG. 4: [color online] Differential asymmetry A(Ay t f) at par- 
ton level for different MC samples. See Ref. [4l|] for the details 
on axigluon models. The observed A{Ay t f) is fitted with the 
functional form of Eq. (0. 


This procedure preserves the even part of the distribu¬ 
tion of Ay t f. It also preserves the original shape of the 
differential asymmetry, but changes its magnitude to the 
desired value. 


3. Calibration 


Starting from the nominal MC@NLO ti sample, we pro¬ 
duce test samples using the product of the weights de¬ 
fined in sections I VIA 11 and IVI A 21 We use a grid of 
values for polarizations of nP = (—0.2,—0.1, 0,0.1,0.2) 
and asymmetries of A tt = (—0.1,0,0.1,0.15,0.2,0.25) to 
obtain 30 samples in addition to the unweighted nom¬ 
inal sample. We apply the method of ME reconstruc¬ 
tion to each of the 31 fully simulated samples and ex¬ 
tract a raw measurement (A“ w ,reP raw ) associated with 
a given parton-level (A u ,kP). A fit to the obtained set 
of points in the space (A u , kP, A‘ aw , /«P raw ) determines 
two affine functions that relate the reconstructed quan¬ 
tities to the true quantities: A( aw = fi(A u ,nP) and 
«P raw = f 2 (A tt ,nP). The affine functions fit the 31 
points well, with residuals < 0.1%. We rewrite the affine 
relations using a matrix equation: 

{«t)= C i«r) +0 ’ (12 > 

where C is a 2 x 2 calibration matrix and O is a vector 
of offset terms. The values of the matrix C and O are 
reported in Table 0 for the the different dilepton chan¬ 
nels. To determine the statistical uncertainties on the 
calibration parameters, we use an ensemble method. We 
split the MC@NLO samples into 100 independent ensem¬ 
bles and then repeat the calibration procedure for each 
of them. 
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TABLE V: Calibration parameters and their statistical uncertainties for the different channels. 


Channel 

Calibration matrix C 

Offset O 

e/i | 

( 0.617+0.008 
v 0.346+0.008 

0.148 + 0.002 \ 
0.560 + 0.003 ) 

( 0.011 + 0.002 \ 

{ -0.009 + 0.003 ) 

ee | 

( 0.599+ 0.006 
^ 0.315+0.007 

0.135 + 0.003 \ 
0.544 + 0.005 ) 

( 0.007+0.003 \ 

i, -0.003 + 0.005 ) 

y/i | 

( 0.639+0.007 
( 0.460+0.008 

0.189 + 0.004 \ 
0.649 + 0.006 ) 

-0.006 + 0.007 ) 

Dilepton | 

( 0.617+0.008 

V 0.359+ 0.006 

0.153 + 0.002 \ 
0.572 + 0.002 ) 

( 0.010 + 0.002 \ 
i, -0.007+0.003 ) 


B. Measurement of A u and kP after calibration 

The calibration relation of Eq. m. is inverted to re¬ 
trieve the true partonic asymmetry A u and the true po¬ 
larization kP from the reconstructed Al* and «P r - w . 
We obtain a measurement of (A tt ,nP) reported in Ta¬ 
ble IVII for each dileptonic channel using the calibration 
coefficients from Tablc[V]and the raw measurements from 
Tables [TT] and Mil 

Two ALPGEN+PYTHIA tt samples generated at dif¬ 
ferent rrit are used to estimate the dependence of the 
measurement on m. t . Considering a top mass of rrit = 
173.18 ± 0.94 GeV [4(J as reference, the dilepton results 
reported in Tabic IVTl have to be corrected by —0.02% 
and 0.15% for A u and kP, respectively. The corrected 
combined dilepton results are 


= (15.0 + 6.4 (stat))% 

(13) 

kP = (7.2 ± 10.5 (stat))%. 

(14) 


VII. SYSTEMATIC UNCERTAINTIES 

We consider three categories of uncertainties. Un¬ 
certainties affecting the signal are obtained by deriving 
calibration coefficients from alternate signal models and 
propagating them to the final results. Uncertainties af¬ 
fecting the background have an impact on the raw mea¬ 
surements, A* aw and K-P raw , as these observables are ob¬ 
tained after subtracting the background. They are prop¬ 
agated to the final measurement by applying the nominal 
calibration correction to the modified At* and 
The third category consists of the uncertainties on the 
calibration method. Since the measurement is performed 
after background subtraction, the calibration is indepen¬ 
dent of the normalization of the tt simulation, and there 
is no systematic uncertainty due to signal normalization. 
The uncertainties on A u and kP due to the different 
sources are summarized in Table IVIII together with the 
correlations. 


A. Uncertainties on signal 

Several sources of systematic uncertainties due to the 
detector and reconstruction model affect the jets and thus 
the signal kinematics. We consider uncertainties on the 
jet energy scale, flavor-dependent jet response, and jet 
energy resolution [32:] . We also take into account uncer¬ 
tainties associated with b tagging and vertexing [34j. 

To estimate the impact of higher order correc¬ 
tion, we compare the calibration obtained with 
MC@NLO+HERWiG to the calibration obtained with 
ALPGEN+HERWIG. To propagate uncertainty on the 
simulation of initial state and final state radiations 
(ISR/FSR), the amount of radiation is varied by 
scaling the ktfac parameter either by a factor of 
1.5 or 1/1.5 in an ALPGEN+PYTHIA simulation of tt 
events [501 ]. The hadronization and parton-shower model 
uncertainty is derived from the difference between the 
PYTHIA and HERWIG generators, estimated by comparing 
ALPGEN+HERWIG to ALPGEN+PYTHIA tt samples. The 
different models for parton showers used by various MC 
generators yield different amounts of ISR between for¬ 
ward and backward events j57j, [58|. The uncertainty on 
the ISR model is defined as 50% of the difference be¬ 
tween the nominal results and the results derived from 
a MC@NLO simulation in which the dependence of the 
forward-backward asymmetry on the pr of the tt s yst em 
is removed. The uncertainty of 0.94 GeV on m* [40| is 
propagated to the final result using two ALPGEN+PYTHIA 
samples generated with different m t values. We deter¬ 
mine PDF uncertainties by vary ing the 20 parameters 
describing the CTEQ6M1 PDF [35( within their uncer¬ 
tainties. 


B. Uncertainties on background 

The uncertainty on the background level is obtained by 
varying the instrumental background normalization by 
50% and the overall background normalization by 20%. 
The model of the instrumental background kinematics 
is varied, using the same method as in Ref. [r|. We 
reweight the reconstructed Ay, cos 9 + , and cosd~ dis¬ 
tributions by a factor of 1 + e x Oband; where dband is 
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TABLE VI: Measurements of A u and kP for each dileptonic channel corrected for the calibration (for rrit = 172.5 GeV). The 
statistical correlation between the two measurements arises both from the statistical correlation of the experimental observables 
and the correction for the calibration. 


Channel 

A 44 (%) 

kP (%) 

statistical 
correlation (%) 

e/i 

11.6 ± 7.8(stat) 

12.6 ± 13.0(stat) 

-48 

ee 

26.1 ± 15.2(stat) 

17.5 ± 26.0(stat) 

-58 

pp 

17.8 ± 16.7(stat) 

—22.2 ± 24.6(stat) 

-52 

Dilepton 

15.0± 6.4(stat) 

7.0 ± 10.5(stat) 

-50 


the statistical uncertainty band of the distribution and 
e = ±1 is chosen to be positive for Ay > 0, cos 9 + > 0, 
cosd~ < 0, and negative for Ay < 0, cos 9 + < 0, and 
cos# - > 0. 


C. Uncertainties on calibration 

We also consider sources of uncertainties affecting the 
calibration procedure. The statistical uncertainty on the 
calibration parameters and their correlations are propa¬ 
gated to the final measurements. The uncertainties are 
0.60% for A tt and 0.61% for kP. The correlation is 
-39%. 

To estimate a systematic uncertainty due to the choice 
of A y t t calibration procedure, we reweight the MC@NLO 
sample to reproduce the shape of the differential asym¬ 
metries of each different BSM and SM model consid¬ 
ered. Each of the resulting samples serves as a seed for a 
new calibration procedure as described in section fVI A 21 
The maximum variation in the A u measurement ob¬ 
tained with these new calibrations is taken as system¬ 
atic uncertainty. It is obtained using the shape from the 
ALPGEN + PYTHIA sample and amounts to 1.3%. The 
impact of these tests is negligible for nP since only the 
A y t i distribution is modified. 

We also perform a closure test using the five differ¬ 
ent BSM models described in section IIV Bl For each 
of the considered BSM models we create test samples 
by reweighting the A y t i and cos distributions, in the 
same way as described in section IVI Al for MC@NLO sam¬ 
ples. The samples cover a range of values of A u and kP 
centered around the data measurement within ±1 statis¬ 
tical standard deviations. These samples are treated as 
pseudo-data: We compute the differences between what 
would be measured using the nominal calibration and 
the true A tt and kP of each sample. The maximum A u 
bias is found for the axigluon m200L sample [41] and 
corresponds to a shift of (AA 44 ,AkP) = (—2.9%, 2.3%) 
obtained for (A 44 ,/cP) « (19%, 9%). The maximum 
kP bias is found for the axigluon m200A sample |4lj| 
and corresponds to (AA 44 ,AkP) = (—1.5%, 2.6%) 

for (A 44 ,kP) ss (10%,0%). These two doublets in 
(AA 44 ,AkP) are taken as uncorrelated systematic un¬ 
certainties. In each of these doublets, the uncertainty on 
A u and kP are taken as —100% correlated. 


VIII. RESULTS 

The measurements and the uncertainties discussed in 
the previous sections are summarized by 

A™ = (15.0 ± 6.4 (stat) ± 4.9 (syst))%, 
kP = (7.2 ± 10.5 (stat) ± 4.2 (syst))%, (15) 

with a correlation of —56% between the measurements. 
The results are presented in Fig. [5] The NLO SM predic¬ 
tion for A tl is A 44 = (9.5±0.7)% @, while the SM polar¬ 
ization is expected to be small, kP = (—0.19±0.05)% Q. 
Our measurement is consistent with the SM prediction 
within the 68% confidence level region. In Fig. [5] we over¬ 
lay the expected values for the different axigluon models 
of Ref. [44]]. As the models are generated with the LO 
MADGRAPH generator, we add an asymmetry of 9.5% 
arising from the pure SM contributions that is not ac¬ 
counted for by MADGRAPH. The approximation of just 
adding the MADGRAPH LO asymmetry to the SM asym¬ 
metry is estimated to be valid at the ~ 3% level. 
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FIG. 5: [color online] Two dimensional visualization of the 
A 44 and kP measurements and comparison with benchmark 
axigluon models p~3| . 

We interpret the measurements as a test of the SM, 
separately assuming the SM forward-backward asymme¬ 
try of A 44 = (9.5 ± 0.7)% and the SM polarization of 
kP = (—0.19 ± 0.05)%. As we assume the SM, we do 























14 


TABLE VII: Summary of systematic and statistical uncertainties. 


Source of uncertainty 

Uncertainty on 

A u (%) Uncertainty on 

kP (%) Correlation (%) 

Detector modeling 




Jet energy scale 

0.13 

0.50 

-100 

Jet energy resolution 

0.03 

0.06 

100 

Flavor-dependent jet response 

0.02 

0.06 

-100 

b tagging 

0.14 

0.43 

- 94 

Signal modeling 




ISR/FSR 

0.16 

0.41 

-100 

Forward/backward ISR 

0.10 

0.07 

-100 

Hadronization and showering 

3.28 

1.94 

-100 

Higher order correction 

0.02 

0.71 

-100 

PDF 

0.12 

0.30 

- 98 

Top quark mass 

0.03 

0.21 

-100 

Background model 




Instrumental background shape 

0.16 

0.53 

100 

Instrumental background normalization 

0.29 

0.01 

100 

Background normalization 

0.44 

0.18 

100 

Calibration 




A y t t model 

1.28 

0.11 

100 

MC statistics 

0.60 

0.61 

- 39 

Model dependence 




Maximum A u variation 

2.91 

2.35 

-100 

Maximum kP variation 

1.49 

2.58 

-100 

Statistical uncertainty 

6.40 

10.53 

- 50 

Total systematic without model dependence 

3.62 

2.40 

- 74 

Total systematic 

4.88 

4.24 

- 83 

Total 

8.05 

11.35 

- 56 


not consider the uncertainty from the dependence on the 
physics model. The constraint on A tt is applied to the 
two-dimensional result of Eq. (fl5l) to obtain the polar¬ 
ization 

kP = (11.3 ± 9.1 (stat) ± 1.9 (syst))%. (16) 

This result is consistent with the SM expectation at the 
1.2 standard deviation level. Applying the constraint on 
kP we obtain an asymmetry of 

A tl = (17.5 ± 5.6 (stat) ± 3.1 (syst))%, (17) 

which is consistent with the SM expectation at the 1.3 
standard deviation level. 

In a previous publication, the DO Collaboration has 
measured the forward-backward asymmetry in the lep- 
ton+jets channel Q 

A tl = (10.6 ± 2.8 (stat) ± 1.3 (syst))% = (10.6 ± 3.0)%. 

(18) 

This lepton+jets measurement was performed in the con¬ 
text of a test of the SM, as no study of the depen¬ 
dence with respect to the possible polarization was per¬ 
formed. Therefore, it should be compared with the re¬ 
sult of Eq. m- We classify the systematic uncertain¬ 
ties of both measurements by their sources and consider 


them as being either completely correlated, e.g ., the b- 
tagging uncertainty, or completely uncorrelated, e.g ., the 
background modeling. Even if some sources of uncertain¬ 
ties are correlated between both channels, the dominant 
sources are not, so that the final overall uncertainties are 
only 7% correlated. The two measurements are consis¬ 
tent with a probability of 30% given by a x 2 test- We 
combine the lepton+jets and dilepton measurements .us¬ 
ing the best linear unbiased estimate (BLUE) [53, |6G |. 
The combination is a weighted average of the input mea¬ 
surements, with the dilepton measurement given a weight 
of 0.17 and the lepton+jets measurement a weight of 0.83. 
The combined result is 

A tl = 11.8 + 2.5 (stat) ±1.3 (syst))%. (19) 


IX. SUMMARY 


We have presented a simultaneous measurement of the 
forward-backward asymmetry of tt production and the 
top quark spin polarization in the beam basis in dilepton 
final states, using 9.7 fb -1 of proton-antiproton collisions 
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at y/s = 1.96 TeV with the DO detector. The results are: 

A tl = (15.0 ±8.0)%, (20) 

kP = (7.2 ±11.3)%, 

with a correlation of —56% between the measurements. 
They are consistent with the SM expectations within the 
68% confidence level region. 

Interpreted as a test of the SM and assuming the SM 
forward-backward asymmetry, these results yield a mea¬ 
surement of the top polarization of 

kP = (11.3 ±9.3)%. (21) 

Assuming the SM polarization, we obtain a forward- 
backward asymmetry of 

A“ = (17.5 ±6.3)%. (22) 

This asymmetry is combined with the measurement of 
the asymmetry in lepton±jets final states yielding a com¬ 
bined asymmetry of 

A a = (11.8 ±2.8)%. (23) 

All of these results are consistent with the SM expec¬ 
tations within uncertainties. 
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